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Imaging with a layered superlens is a spatial filtering operation characterized by the point spread function 
(PSF). We show that in the same optical system the image of a narrow sub- wavelength Gaussian incident 
field may be surprisingly dissimilar to the PSF, and the width of PSF is not a straightforward measure of 
resolution. FWHM or std. dev. of PSF give ambiguous information about the actual resolution, and imaging 
of objects smaller than the FWHM of PSF is possible. A multiscale analysis of imaging gives good insight into 
the peculiar scale-dependent properties of sub-wavelength imaging. © 2010 Optical Society of America 
OCIS codes: 160.4236 , 050.0050, 110.0110 , 260.0260 , 310.6628, 100.6640 
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Fig. 1. Intensity transmission and reflection coefficients 
(T, R) of the multilayer as a function of the filling fac- 
tor, plotted over the transfer function (in vertical cross- 
sections of the color map). Phase isolines are distanced 
by 7r/4. 
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Fig. 2. Point spread function (in vertical cross-sections) 
in the function of the filling factor (top) compared with 
the response of the multilayer to a narrow Gaussian sig- 
nal with sub- wavelength width (a(Hy) = 0.2 A) (bot- 
tom). Phase isolines are distanced by tt/2. 



Since the seminal paper by Pendry [1] and the first 
experimental demonstration of sub- wavelength imaging 
through a 40nra silver slab [2,3], sub- wavelength imag- 
ing at visible wavelengths has been investigated in much 
thicker low- loss layered silver-dielectric structures [4-10]. 
A variety of physical models may be applied to explain 
the mechanism of transmission: 1. the effective medium 
anisotropic approximation of the sub-wavelength multi- 
layer [4] combined with the Fabry-Perot resonant con- 
dition tuned to be independent of the angle of inci- 
dence [8,9]; 2. negative refraction at the boundary of 
silver slab [1] multiple interfaces of the layered structure 
resulting in diffract ion- free propagation [6]; 3. resonant 
tunneling through the bandgap material formed by the 
periodic metal-dielectric multilayer [5]. Enhancement of 



evanescent fields needed for sub- wavelength imaging may 
be also explained in various ways: 1. as the result of col- 
lective coupling between plasmon modes at subsequent 
metallic layers [11, 12] - if we look to the internal field 
distributions; 2. as self-collimation [13] - if we look to the 
band diagrams of the multilayer; 3. as the result of large 
effective permittivity e± « oo - when the homogenized 
anisotropic model of the system is valid. 

Here, we calculate numerically the transfer function 
(TF) of the multilayer system. Let us emphasize that 
the structures analyzed are inaccurately characterized 
with the homogenization model which underestimates 
the value of transmission. 

Imaging through a layered superlens consisting of uni- 
form and isotropic materials is a linear and shift invari- 
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Fig. 3. TF (top) and PSF (bottom) of the two multilayers 
(solid line- amplitude/intensity; dashed line -phase). 
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Fig. 4. Transmission of a sub- wavelength Gaussian inci- 
dent field of various width through two multilayers (ver- 
tical cross-sections include the normalized field H 2 ). 

ant (LSI) spatial filtering operation. Furthermore, pla- 
nar imaging is a scalar LSI system. The theory of LSI 
systems, widely utilized in the diffraction theory and in 
Fourier Optics is therefore readily applicable. However, 
for systems with sub- wavelength resolution, point spread 
function (PSF) engineering leads to surprising peculiari- 
ties. In this paper we show that the width of PSF is not 
a straightforward measure of resolution. 

For the TM polarization, the magnetic field H y (x,z) 
may be represented with its spatial spectrum H y (k xi z) 

/+oo 
H y (k x , z)exp{ik x x)dk x , (1) 
-co 

where, at least for lossless materials, the spatial spectrum 
is clearly separated into the propagating part k 2 < k^e 
and evanescent part k 2 > k^e. 



The TF is the ratio of the spatial spectra of the out- 
put and incident fields and corresponds to the amplitude 
transmission coefficient of the multilayer 

H y (k x , z = L)= TF(k x ) ■ Hl nc (k x , z = 0). (2) 

Due to reflections, the incident field H y nc (k x ,z = 0) 
differs from the total field H y (k xi z = 0). 

PSF of an imaging system can be often straightfor- 
wardly interpreted, and provides clear information about 
the resolution, loss or enhancement of contrast, as well 
as the characteristics of image distortions. For instance, 
the resolution may be usually linked to the width of PSF. 
When the input signal and PSF are positive functions 
with a limited region of support (region with non-zero 
values), the regions of support of convolved functions 
simply add together, contributing to the broadening of 
the filtered signal. This can be expressed more conve- 
niently using the following relation on the Lq norms, 

\\H y (x) * PSF(x)\\ = \\H y (x)\\o + \\PSF(x)\\o, (3) 

where ||/(x)|| = lim p ^ + f\f(x)\ p dx. 

On the other hand, for simple Gaussian PSF and in- 
put, the output has the width (variance) equal to the 
sum of variances of PSF and input, 

exp(— x 2 /cr 2 ) *exp(— x 2 /cr 2 ) oc exp(— x 2 j(<j\ + cr|)). (4) 

Therefore once again the width of PSF has a clear link 
to the resolution of the imaging system. However, formu- 
las (4) and (3) are often invalid for diffractive systems 
with complex PSF. 

From now on, we focus on a SrTiOs — Ag multi- 
layer with N = 20 periods, and the total thickness 
L = 1.15/im. The elementary cell consists of an Ag 
layer symmetrically coated with SrTiOs. Strontium Ti- 
tanate is an isotropic material with a high refractive in- 
dex n = 2.674+0.027z at the wavelength A = 430nra [14]. 
The refractive index of silver at the same wavelength is 
equal to ua q — 0.04 + 2.46z [15]. Fig. 1 shows how the 
TF of the multilayer depends on the silver-filling frac- 
tion, and the corresponding PSF is shown in the up- 
per part of Fig 2. The evanescent part of the TF has 
a large magnitude, which is the necessary condition for 
sub-wavelength imaging. The shape of TF is generally 
regular with the exception of the phase discontinuity in 
the vicinity of k x /ko = 1, as well as the strong phase 
modulation below dA g /A < 0.35 which suppresses the 
super-resolving properties of the PSF in that range. The 
phase step at k x /ko = 1 in the TF influences the shape of 
the corresponding PSF which, with the increase of filling- 
factor, evolves from a narrow sub- wavelength maximum 
to a shape dominated by the side-lobe. The response 
to a narrow sub- wavelength Gaussian signal is entirely 
different from the PSF (Fig2, bottom). PSF does not 
resemble a Gaussian function and its width measured 
with FWHM is different from the doubled standard de- 
viation. The off-axis background of PSF results in the 
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Fig. 5. Transmission of light through multilayer I (left) 
and multilayer II (right) (Instantaneous distribution of 
\H y \ obtained with FDTD). The incident field is a point 
source (top) or a narrow Gaussian object with cr(Hy) = 
0.2A (bottom) 



high value of std. dev., and probably FWHM is a more 
meaningful measure of resolution of the system. More- 
over, the broadening of the optical signal can not be 
expressed with formulas (4) or (3). The exception is the 
range of filling fraction in between 0.35 and 0.45, where 
the PSF resembles a Gaussian function and the broad- 
ening follows a simple intuitive behavior. 

More in general, the width of response may even show 
an anomalous (decreasing) dependence on the size of the 
sub-wavelength Gaussian incident signal. It is striking 
how dissimilar are the PSF and the response to a narrow 
Gaussian signal around cUp/A ^ 0.5 (Fig 2). The expla- 
nation is nevertheless not difficult, as the bandwidth of 
the TF in use depends (inversely) on the width of the in- 
cident Gaussian function. The opposite phase of TF for 
propagating and evanescent waves is the source of the 
side- lobe of PSF. Partial removal of the central maxi- 
mum in PSF (equal to the mean value of TF) occurs 
only when the contribution from evanescent and prop- 
agating harmonics to the mean value compensate each 
other. Broader Gaussian incident fields limit the band- 
width in use, and suppress this sensitive condition. 

In Figs. 3, 4 we show the TF, and PSF of two se- 
lected multilayers with the filling fractions equal to 0.53, 
and 0.37, as well as their response to a sub-wavelength 
Gaussian field distribution with FWHM < 1.6A. These 
two multilayers represent the situation of a regular 
nearly Gaussian PSF and a side-lobe dominated PSF, 
respectively. Both multilayers allow for imaging of sub- 
wavelength details, however their responses scale differ- 
ently with the size of sub- wavelength object. We at- 
tribute their different behavior to the value of the phase 
shift between the propagating and evanescent part of 
the TF (fig. 3). We have recently analyzed an analogous 
situation [10], however resulting from the different mod- 
ulation depth of TF. 

Finally, let us demonstrate the same example with 
FDTD simulation obtained with the open source Meep 
package [16]. In Fig. 5 we show how a narrow Gaussian 



beam, as well as a beam originating from a point source 
propagate through the two discussed multilayers. The 
simulations only confirm the behavior described in Fig. 4 
but are a good illustration to the peculiarities encoun- 
tered in sub- wavelength imaging. While one of the mul- 
tilayers allows for approximately diffraction-free propa- 
gation independently of the size of the source, the other 
behaves in the same way for broader sources only, and 
shows strong diffraction when the shape of the source 
approaches a (^-function. 

In conclusion, we have studied the transmission of sub- 
wavelength incident Gaussian field through a thick (L ~ 
nX) silver-dielectric superlens. We have demonstrated 
that the response to narrow sub-wavelength Gaussian 
signal may be surprisingly different from the PSF of the 
system. Multiscale analysis provides the means to dis- 
tinguish between diffraction-free propagation for various 
ranges of object sizes. 

We acknowledge support from the Polish projects 
N N202 033237, N R15 0018 06 and the framework of 
COST actions MP0702, MP0803. 
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